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We report on the magnetic penetration depth, A, in a type II superconductor NbB2+i: 
determined by muon spin rotation method. We show in the sample with a; = 0.1 that A at 
2.0 K is independent of an apphed magnetic field. This suggests that the superconducting 
order parameter in NbB2+a; is isotropic, as expected for conventional BCS superconductors. 
Meanwhile, the superconducting carrier density (oc A~^) exhibits an interesting tendency 
of increase with increasing (where varies with x). Possible origin of such behavior is 
discussed in comparison with the case of exotic superconductors. 
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The recent revelation of superconductivity in MgB2 with high transition temperature (Tc = 
39K)^ has revived much interest in hexagonal diborides which were once subjected to active 
studies from 1950's through 70's.^"^ The discovery of MgB2 prompted various theoretical 
and experimental investigations to clarify the origin of the high Tc seemingly exceeding the 
so-called 'BCS-barrier' (~ 30 K). These compounds have a simple layered structure along 
c-axis, consisting of an alternative stack of metal ion layers forming a triangular lattice and 
boron atom layers in a honeycomb network. In MgB2, both theoretical and experimental 
studies suggested that the strong electron-phonon interaction in accordance with the light 
mass of the organized atoms and two-dimensional character of boron atom network are closely 
related to the relatively high-Tc superconductivity. Moreover, it turned out that some of the 
unconventional properties in the superconducting phase can be attributed to the multiband 
superconductivity, where the presence of two energy gaps corresponding to p-,^ and po- orbits 
of boron atoms have been confirmed by experimental^ and theoretical studies.^ 

Among various diborides, NbB2 is intriguing as it exhibits strong sensitivity of supercon- 
ducting transition temperature to the stoichiometric imperfection; Cooper et al. reported Tc 
= 3.87 K for boron-rich NbB2,^ while Leyarovska et al. reported much lower Tc ( = 0.62 K) 
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for a stoichiometric NbB2.^ Schirbcr et al. showed = 9.4 K for single-crystaUine sample of 
NbBj; (x ~ 2).*^ Recently, Yamamoto et al. reported that changes between 2 K and 9.2 K 
in the samples prepared under high pressure, probably because of the shift in the composition 
ratio between Nb and boron. It is noticeable that Escamilla et al. found a maximum = 9.8 
K,^ which is slightly higher than that of elemental Nb (=9.23 K). Thus, it might be possible to 
obtain a guide for ever higher in hexagonal diboridcs by studying the relationship between 
Tc and various parameters including chemical composition and crystal structure in NbB2-f-a;* 

In this paper, we report on the magnetic penetration depth. A, in the mixed state of 
NbB2+x studied by muon spin rotation and relaxation method (^uSR). The magnetic penetra- 
tion depth, which reflects superconducting carrier density, rig (oc A~^), is determined micro- 
scopically by yuSR signals due to the inhomogeneity of magnetic field distribution in the flux 
line lattice (FLL). We show that A is independent of the applied field, strongly suggesting that 
NbB2+a; has an isotropic superconducting order parameter. More interestingly, a quasi-linear 
relationship has been observed between the muon spin relaxation rate a (T = K, which is 
proportional to rig) and Tc as it varies with x. This result suggests that the change of in 
NbB2+x depends on the superconducting carrier density. 

Polycrystalline samples of NbB2+x (x = 0, 0.03, 0.1 where x is a nominal value) used in 
this experiment were synthesized from Nb powder (99.96 %) and boron crystalline powder 
(99 %). The pellets of these mixtures were sealed by quartz tube at vacuum pressure (< 
3.0 X 10^^ Torr) and heated at 1000 °C for 30 hours. These samples were confirmed to be 
single phase, and to decrease a-axis and incrcse c-axis with increasing ,x by the X-ray powder 
diffraction pattern using the conventional X-ray powder spectrometer (RAD-C: RIGAKU). 
The magnetic susceptibility and electrical resistivity were measured with the SQUID magne- 
tometer (MPMS2; Quantum Design Co., Ltd.) and the PPMS system (Quantum Design Co., 
Ltd.). 

Figure 1 shows the temperature dependence of magnetic susceptibility upon field cooling 
at 10 Oe, where one can observe significant reduction of the susceptibility at a well-defined 
temperature, Tc, for respective compositions. This is obviously due to the Meissener dia- 
magnetism associated with the superconducting transition. It is also clear in Fig. 1 that 
increases with increasing x. The superconducting volume fraction at 1.8 K is estimated to 
be approximately 44, 34 and 29 % for NbB2, NbB2.o3 and NbB2.i. The corresponding upper 
critical field (i7c2(0)) are 6.4, 7.6 and 10.5 kOe respectively, which were estimated from the 
results of electrical resistivity measurements under an applied magnetic field. 

The /iSR experiment was performed on the M15 and M20 beamlines at TRIUMF which 
provides a muon beam with the momentum of 29 MeV/c. The polycrystalline samples were 
mounted on sample holders having a dimension of 7 X 7 mm^ (M15) or 25 x 25 mm^ (M20). 
The yuSR spectra were obtained at H = 1 kOe and at T = 2 K to map out the temperature 
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and magnetic field dependence, respectively. Care was taken to cool down the sample in the 
field-cooled condition at each applied field to avoid the effect of flux pinning. 

Since the muons stop randomly on the FLL, the muon spin precession signal P{t) provides 
a random sampling of the internal field distribution B{f), 

/oo 
n{B)e^v{ilnBt)dB, (1) 
-oo 

n{B) = {6{B{f)-B))r, (2) 

where 7^ is the muon gyromagnetic ratio (= 2ttx 13.553 MHz/kOe), and n{B) is the spectral 
density for muon precession determined by the local field distribution. These equations indicate 
that the real amplitude of the Fourier transformed muon precession signal corresponds to 
the local field distribution n[B). The London penetration depth in the FLL state is related 
to the second moment ((Ai?)^) = {{B{r) — H)'^) of the field distribution reflected in the 
^SR line shape. In poly crystalline samples, a Gaussian distribution of local fields is a good 
approximation, 

P{t) ~ exp(-a2iV2)exp(n^m), (3) 

a = 7M\/((ABp). (4) 

For the ideal triangular FLL with isotropic effective carrier mass m*, A is given by the rela- 
tion,io-i2 

a [fis-^] = 4.83 X 10^(1 - h)[l + 3.9(1 - hfj^/^X'^ [nm], (5) 
where h = H/Hc2, and A is related to the superconducting carrier density, 

A' = 3^. (6) 

47rnse"' 

indicating that A is enhanced upon the reduction of Us due to the quasiparticle excitations. 

Figure 2 shows typical muon spin precession signals in NbB2.i at various temperatures 
under a transverse field of 1 kOe, where one can observe the change in the lineshape due to FLL 
formation as the temperature passes through Tp. The damping of precession amplitude above 
Tc is mainly due to static random local fields from ^^Nb and ^'^'^^B nuclear magnetic moments. 
Upon cooling below Tc, the spectrum exhibits enhanced depolarization, which indicates that 
magnetic field distribution becomes inhomogeneous due to the formation of magnetic vortices. 
In analyzing these spectra, we adopted the following fitting function, 

where the index j denotes the components of normal {j = 1) and superconducting {j = 
2) domains, A is the total positron decay asymmetry, Aj is the partial asymmetry, aj is 
the muon spin relaxation rate, Bj is the central frequency and is the initial phase for 
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respective components. We define ai to reflect the spin relaxation in the normal state where 
it is dominated by the nuclear magnetic moments, dn, and (T2 in the superconducting state 
where o"! = 0"„ + with being due to the formation of FLL. The temperature dependence 
of cJv in NbB2.i a,t H = 1 kOe is shown in Fig. 3. In the FLL state, increases with decreasing 
temperature below rc(l kOe) ~ 5 K. According to the empirical two-fluid model, which is 
supposed to be a good approximation for the BCS theory, the following relation is expected, 

A(T) = A(0)^=i=, (8) 

which leads to 

^7.(r)=a.(0)(l-T^), (9) 

where r = T/Tc(l kOe). As clearly seen in Fig. 3, the agreement between the above relation 
and our data is far from satisfactory. A better agreement is attained when we perform the 
fitting analysis by a similar formula, 

(7v(r) = (7v(0)(l-r2). (10) 

It yields Tc = 5.3(1) K when Tc is a free parameter, which is in good agreement with the value 
estimated from the electrical resistivity measurement under the same field 5.0 K). We 
also performed fitting analysis by a weak-coupling BCS model (w-BCS),^^ which turned out 
to reproduce the experimental data very well. Thus, these analyses indicate the difficulty to 
assess the structure of order parameter in NbB2+a; based solely on the temperature dependence 
of (Tv with the data above 2 K. 

The situation is much improved by looking into the magnetic field dependence of a^,. As 
shown in Fig. 4(a) in NbB2.i decreases with increasing external field, where the solid line 
represents the fitting result by eq.(5) with Hc2 = 6.7 kOe determined by electrical resistivity 
measurements. The corresponding field dependence of A is shown in Fig. 4(b), where A is 
mostly independent of applied magnetic field. For a quantitative evaluation on the strength 
of the pair-breaking effect by external field, we performed a fitting analysis by the following 
simple linear relation, 

A(/i) = A(0)[l + 77/i], (11) 

where 77 is a dimension-less parameter which represents the strength of pair breaking effect, and 
h = iJ/i?c2(2 K) with i?c2(2 K) = 6.7 kOe. Prom the fitting analysis for the data by eq.(ll), we 
obtain 77 = 0.02(2) with A(0) = 233(2) nm, where the result is shown as a solid line in Fig. 4(b). 
It has been established experimentally that the parameter 77 represents the degree of anisotropy 
in the superconducting order parameter,^^ since the pair-breaking mechanism due to the quasi- 
classical Doppler shift has come to wide recognition.^^ For example, in Y(Nio.8Pto.2)B2C,^^ 
Cd2Re207^^ and ¥38! {h < 0.5),^^ which behave as conventional BCS superconductors with 
isotropic s-wave pairing, 77 is reported to be nearly zero. On the other hand, r/ = 5 - 6.6 
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in a high-Tc cuprate YBa2Cu307_5 which has a d-wavc pairing. In the case of YNi2B2C 
which turned out to be an anisotropic s-wave superconductor, r/ ~ 1 is reported. ^'^ From 
the comparison of these results with our result, it is strongly suggested that NbB2.i has an 
isotropic s-wave order parameter. 

We made additional measurements for NbB2+x with a; = and 0.03 to clarify the relation 
between the bulk and microscopic parameters such as A. The relation between and 
(Tv(0), which is extrapolated from the temperature dependence of a^{T), is shown in Fig. 5. 
By the comparison between eqs.(5) and (6), the following relation is derived; 

CTv oc -2 oc — . (12) 

It is inferred from Fig. 5 that Tc increases with increasing Us/m* with a quasi- linear relation 
(Tc oc ns/m*) . This strongly suggests that the change of Tc is associated with that of the 
superconducting carrier density. It is interesting to note that this behavior has a certain simi- 
larity with the case of exotic (e.g. high-Tc cuprate, organic, heavy-fermion) superconductors.^^ 
According to the MEM/Rietveld analysis on the synchrotron radiation X-ray measurements, 
this corresponds to the change in the number of electrons at Nb- and B-sites.^^ It is generally 
expected from the BCS theory that 

T, oc exp(-l/iV(0)F), (13) 

where A^(0) is the density of state at the Fermi surface and V is the electron-phonon coupling 
constant. Assuming that V is independent of the slight change in the chemical composition, our 
result suggests that A^(0) is strongly dependent of the boron content x in NbB2+a;. This might 
suggest that a similar effect may have to be considered in understanding the Ucmura-plot 
for exotic superconductors.^^ However, we note that the possibility of x-dependent electron- 
phonon coupling cannot be ruled out solely based on the present result, considering the pos- 
sible modulation of lattice structure by non-stoichiometric boron contents. 

In summary, we have investigated the magnetic penetration depth in the FLL state of 
NbB2+a; using /xSR method. The field dependence of A strongly suggests that the supercon- 
ducting order parameter in NbB2+2: is isotropic and thereby described by the BCS s-wave 
pairing. The composition dependence of Tc in this compound suggests that the superconduct- 
ing carrier density is a limiting factor in determining Tc as inferred from the relationship 
between Tc and a^. 
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Fig. 1. Temperature dependence of the magnetic susceptibility at 10 Oe measured under field-cooling 
condition. 
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Fig. 2. Typical muon spin precession signals in NbB2.i {T^ = 5.0 K) under a transverse field of 1 kOe 
at temperatures of (a) 10 K (above T^), (b) 4 K and (c) 2 K (below T^). 
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Fig. 3. Temperature dependence of the muon spin relaxation rate ctv due to flux line lattice formation 
at iJ = 1 kOe. The curves are results of fitting by a relation (Tv(T') = a'(0)[l — (T/Tc)''] and weak- 
coupling BCS model (dot-dashed curve). The solid curve is for /? =2 and the dotted curve is for 
/3 = 4. 
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Fig. 4. Magnetic field dependence of (a) the muon spin relaxation rate CTv and (b) the magnetic 
penetration depth A in NbB2.i at 2 K. The curves in (a) and (b) are results of fitting by eqs.(5) 
with A proportional to A(0)[1 + rih\ and (11). 
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Fig. 5. The correlation between Tc at H = 1 kOe and av{T = K) in NbB2+x, where the latter is 
estimated from the temperature dependence of ay. 
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